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Expoirimental Study of the Thermal-Acoustic Efficiency 
in a Long Turbulent Diffueion-Plame Burner 

1 . INTRODUCTION 

This report describes a two-year study of noise production in 
a long tubular burner. The research was motivated by an interest 
in understanding and eventually reducing core noise in gas turbine 
engines. The general approach is to employ an acoustic soUrce/pro- 
pagation model to interpret the sound pressure spectrum in the 
acoustic far field of the burner in terms of the source spectrum 
that must have produced it. In the model the sources are assumed 
to be due uniquely to the unsteady component of combustion heat 
release; thus only direct combustion-noise is considered. The 
source spectrum is then the variation v;ith frequency of the 
thermal-acoustic efficiency , defined as the fraction of combustion 
heat release which is converted into acoustic energy at a given 
frequency. 

The thrust of the research was to study the variation of the 
source spectrum with the design and operating parameters of the 
burner. The immediate goal, which was only partially realized in 
the present effort, was to develop a universal model for the spectral 
distribution of the thermal-acoustic efficiency which could then be 
used with the acoustic source/ jvopagation model to predict far-field 
sound pressure spectra for a given burner design. The ultimate goal 
of such research would be the definition of a general acoustic source/ 
propagation model applicab. e to a wide range of gas turbine engines . 
The present research is envisioned as a demonstration of the feasi- 
bility of this idea. 
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2. REVIEW OP PREVIOUS COMBUSTION NOISE RESEARCH 

Although much research has been done on combustion noise ^ con- 
flicting theories still exist concerning the dominant generating 
mechanisms. Many theories of combustion noise have been formulated 

p * 

since Lignthill'S'^ aeroacoustics theory was published in 19S2. 
Strahle^'^, for example, developed three different theories on the 
basis of Lighthill's work. Each of these theories is based on 
fluid mechanics principles and physical reasoning concerning tur- 
bulent flame structures. Strahle argues that the far-field sound 
pressure is directly proportional to the first Eulerian time 
derivative of the chemical reaction rate integrated over the 
reacting volume. Due to the inability to explain the direction- 
ality of the source, Strahle concludes that the directionality 
must be explained by other phenomena such as refraction by tem- 
, perature gradients. Each of the theories developed by strahle^ 
is consistent with the monopole behavior and low frequency nature 
of combustion noise. 

In an extension of Lighthill's theory, Dowling*^ was able to 
explain directionality effects by accounting for variations in the 
mean temperatures and the flow velocity. As anticipated by Strahle, 
the refraction effect, which is included through a Green's function, 
gives the directionality effects missing in references 2 and 3 . 

Gumps t;/^ addresses the generation of combustion noise due to 
entropy disturbances. This type of noise generation is caused by 
hot (or cold) fluid elements passing through the pressure gradients 
of the turbine assembly. Cumpsty concludes that the entropy noise, 
which is a form of indirect combustion noise , dominates the low 
frequency spectrum. 
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The views of Strahle^'^ and Cumpsty® as to the dominant gen- 
erating mechanism in combustion noise are in conflict. The support 
for the direct combustion noise generating mechanism is based mainly 
on isolated combustion systems. In a paper* by Mathews, Rekos, and 
Nagel®, a comparison oif the two main combustion noise generating 
mechanisms in isolated combustors shows that the heat release vari- 
ations caused by the combustion process is the dominant source of 
noise radiated from the engine exhaust. However, Cumpsty® indicates 
that the indirect combustion noise prediction method used, which is 
due to Pickett^, is inadequate for correctly predicting the effects 
of the turbine stages. A more recent paper by Muthukrishnan, 

Strahle, and Neale® shows that results based on isolated combustors 
can be misleading. According to their analysis, combustion noise 
generation is mainly due to direct combustion noise at low exit Mach 
numbers, whereas at high exit Mach numbers (M<0,8) entropy or in- 
direct noise dominates over direct combustion noise. 

Experimental research shows that combustion noise occurs mainly 
at low frequencies and has very little directionality. Motsinger 
and Emmerling® show that for four turboshaft engines ranging in size 
from 250 to 22,800 H.P., the peak of the combustion noise spectrum 
occurs between 400 and 500 Hz. Combustion noise spectra for several 
turbofan and turbojet engines also have peaks between 400 amd 500 Hz. 
The directivity peaks at about 120 deg from the inlet axis for all 
of the engines tested. Methods for predicting the combustor sound 
power level, spectrum, and directivity are also given in reference 9. 

Mathews and Rekos determine that the combustion spectra for 
both can-type and annular geometries were similar. Both experi- 
mental and analytical studies are presented. The peak frequency is 
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found to depend upon the type of burner used. For the burners tested, 
the range of frequencies corresponding to the peak lies between 280 
and 500 Hz. The combustion noise directivity pattern obtained from 
four turbofan engines shows only minor directivity effects. For all 
of the engines tested, the directivity peaks at about 120 deg from 
the inlet axis. 

Smith and Kilham^^ provide fundamental data on sound generated 
by premixed open turbulent flames. Sever U different parameters are 
varied to determine their effect on combustion noise. Examination 
of the combustion-noise spectra indicates that the frequencies 
corresponding to the peaks may be expressed as a constant nondimen- 
sional frequency, or a Strouhal number, defined in terms of the exit 
diameter, the flow and combustion velocities and the frequency of 
the peak. At low velocities, the directional effects are very weak. 
However, as the flow rate increases, the directionality of the com- 
bustion noise source increases. In addition, the angle of maximum 
directivity lies between 40 deg and 80 deg from the burner axis, 
which corresponds to 140 deg to 100 deg from the inlet axis. 

Knott provides data on noise generated by turbulent diffusion 
flames. An annular bi’tner is used with the fuel entering through the 
outer annulus and the oxidizer through the center tube. Results are 
presented for several types of fuels. The directional effects are 
shown to be weak for most of the fuels used. Knott also shows that 
the acoustical energy produced by combustible gas jets is consider- 
ably greater than that emitted by the unburned gas jets. 

Shivashankara , Strahle, and Handley^ ^ report the results of 
further investigations of the noise generated by open turbulent 
flames. Scaling laws for the spectral content, the acoustic power 
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and the thermal-acoustic efficiency are developed for premixed flames. 

An empirical relationship for the acoustic power is developed using 
a regression analysis of the experimental data. Similar results are 
obtained for the three different fuels used. The directionality is 
shown to be weak for all of the premixed flames. The spectra of the 
combustion noise are shown to be predominantly of a low frequency 
nature. The broadband combustion noise has a single peak and a 
gradual amplitude roll-off on either side ?f the peak. The use of 
a Strouhal number to nondimensionalize cornlmstion noise peak fre- 
quencies as suggested by Smith and Kilham^^ is shown to be unaccept- 
able. 

Strahle and Shivashankara^^^ have extended their combustion noise 
studies to include gas turbine combustors. The results of experiments 
conducted to determine the acoustic power and spectra emitted from gas 
turbine combustors exhausting to the atmosphere indicate that the 
acoustic power radiated is independent of fuel/air ratio at low fuel/air 
ratios. However, a transition to a linear dependence on fuel/air ratio 
occurs above a fuel flow for which the primary stabilization region 
becomes fuel rich. The spectra are shown to depend mainly on the 
enclosure acoustics . The smoothed combustion spectrum has a marked 
peak at 400 Hz and weak secondary peaks at 1200 and its harmonics. 

Combustion noise characteristics of a can-type combustor are 
presented by Shivashankara and Crouch^ The effects of several 
different parameters on the spectral content, the directionality, 
the radiated sound power, and the thermal-acoustic efficiency are 
given. A relationship between in-duct and far-field sound levels 
is established. Empirical relationships for the acoustic power 
output and the thermal-acoustic efficiency are also developed, and 
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a generalized combustion noise spectrum is given for several test 
conditions. The generalized spectrum peaks at 315 Hz while the range 
of the peak frequency is between 200 and 600 Hz. The generalized 
directionality peaks at 135 deg from the inlet axis. At low exit 
velocity, the directivity is very weak but increases as the exit 
velocity increases. Individual spectra show resonant peaks due to 
the longitudinal modes of the can-type combustor. The thermal- 
acoustic efficiency of the can-type combustor is two to four orders 
of magnitude higher than the thermal-acoustic efficiency of the open 
turbulent flames in reference 9. 

Ms;thews and Rekos^® develop an improved prediction system for 
direct combustion noise in turbopropulsion systems. Expressions 
for peak frequency and acoustic power output are derived based on 
burner design and performance parameters. Comparison of analytical 
models with experimental data show good correlation. Generalized 
curves of the combustion noise spectrum and directivity show 
similar results to those obtained in references 9, 10, 14, and 15. 

In addition, resonance peaks also appear in the individual com- 
bustion noise spectra. 

Strahle and Muthukrishnan^^ report results of tests of fifteen 
combustors on four different rigs to determine the correlation between 
the experimental and theoretical acoustic power. A theoretical basis 
for the correlation is established using the plane wave theory of 
direct and indirect combustion noise. A comparison of the theoreti- 
cal and experimental regression analysis results is given. A single 
correlation formula is shown to be valid over a wide range of acoustic 
power. It is emphasized that this is a combustor rig correlation for 
equivalent power in a reflection-free situation. The thrust of this 
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study is toward combustor rigs; only the sources directly generated 
by the combustion process or the combustor rig are investigated. In 
other words, only the direct combustion noise and the vorticity corn- 
ponent of the indirect combustion noise are considered; entropy 
noise is not present in this investigation. 

Shivashankara^ ® uses a commercially available auxiliary power 
unit (APU) gas turbine engine and internal-to~far-f ield correlations 
to isolate engine core noise sources. Below 400 Hz, direct com- 
bustion noise is found to dominate the APU far-field noise. Addit- 
ional noise generation between the combustion exit and the turbine 
exit is observed near 375 and 600 Hz . The additional noise gener- 
ation is apparently due to entropy noise or flow noise in the turbine 
inlet torus. Within the secondary duct, mixing of the engine exhaust 
with the cooling air contributes significantly to the far-field APU 
noise. 

Strahle, Muthukrishnan, and Neale investigate experimental 
and analytical separation of hydrodynamic, entropy, and combustion 
noise in gas turbine combustors. Both analytical and experimental 
results indicate that at low Mach numbers , the exterior radiated 
sound is mostly due to direct combustion noise while at high Mach 
numbers, entropy noise dominates the direct combustion noise. 
Analytical results predict that hydrodynamic noise should dominate 
combustion noise below ISO Hz. However, experimental investigations 
show that the hydrodynamic noise is only significant below 50 Hz. 
Attempts at separating combustion and entropy noise between zero 
and 100 Hz are fruitless because of the high correlation between 
them in this frequency range. 
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Milas^® uaes a dynamic system state-space approach to study a 
ducted combustion system. Procedures for calculating pressure auto- 
spectra and cross-spectra given. Comparison of the measured 
spectra with the calculated spectra shows good agreement when 
allowance is made for dispersion and attentuation effects and the 
boundary conditions are properly specified. Resonant peaks due to 
the longitudinal duct modes are observed in the autospectra. 

Strahle^^ presents r brief introduction to the physics of 
combustion noise followed by a short Introduction to some of the 
mathematical and physical descriptions of noise. Comparison of 
experiment and theory for a premixed, gaseous, small jet flame sets 
the stage for discussion of more complex flame types . Although the 
physics of the problem are well understood, a lack of knowledge of 
turbulence details in a turbulent flame makes it impossible to pre- 
dict noise output and spectral content. 

Legendre investigates noise generating mechanisms of a turbu- 
lent gas. He relates noise generation to variations in the circulation 
and diameter of vortex tubes. Density variation in the flow due to 
the combustion process is also considered. An estimate for the sound 
source term in turbulent combustion is presented, as is a corresponding 
wave equation for propagation of sound into the far-fi«ld. 

Ho and Doyle develop an empirical expression for the prediction 
of direct combustion noise by correlating the combustion noise with 
combustion steady-state parameters for a variety of gas turbine 
engines. The spectral shape and far-field directivity of the com- 
bustion noise correlate well with data provided in references 9, 10, 

14, 15/ and 16, The empirical formula for the prediction of direct 
combustion noise is found to correlate within 3-5 dB of the experi- 
mental data for all of the gas turbine engines tested. 
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Muthukrishnan , Strahle, and Neala^** daaeribe experlmants con- 
ducted to determine the cause? of gas turbine combustor noise. 
Results indicate that direct combustion noise and entropy noise 
are the only significant sources of far-field noise. Direct com- 
bustion noise predominates at higher frequencies (f<800 Hz) while 
entropy noise predominates at low frequencies. However, the en- 
tropy noise has only minor effects overall. 

The terminating or radiation impedance at each end of the com- 
bustion system is needed as a boundary condition in order to imple- 
ment any mathematical ttit^del for the standing wav© structure. 

Levine and Schwinger^ ^ develop an approximate solution for 
the radiation of sound from a semi-ir-'^’lnite circular pipe with zero 
wall thickness. The difference between this approximation and the 
experimental results for a thin-v/alled pipe at room temperature is 
minor , 

Ando^® develops an approximate solution for sound radiation 
from a semi-infinite circular pipe having a finite wall thickness. 
For a wall thickness of zero, Ando’,s results are in good agreement 
with the results of Levine and Schwinger^®. Ando's approximation 
has the advantage of being able to vary the wall thickness of the 
pipe so that any pipe, regardless of thickness, can be analyzed. 
Neither the Ando nor the Levine and Schwinger approximation takes 
into account the temperature dependence of the terminating imped- 
ance. Since the present investigation involves implementation of 
the terminating impedance boundary conditions when the exhaust gas 
temperature is elevated siiignificantly above the ambient temperature, 
it is important to understand the effect of this temperature differ- 
ence on impedance. Several attempts have been made to determine the 
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affect of t«?|!poratui:o difference on the terminating impedance; how-* 
ever, investigators have yet to determine the relationship between 
the temperature and the terminating impedance. 

Frickar and Roberts^'^ use the classical impedance tube techniqu® 
to study the acoustic impedance of the open end of a thick**walled 
tube with hot flow. They show the radiation impedance as a function 
of both a wavelength parameter ka and a dimensionless frequency k^a. 
The wave number k in the wavelength parameter depends on both the 
frequency and the temperature of the gas flowing from the pipe, while 
the wave number k^ in the frequency parameter is evaluated at a con- 
stant temperature of 0®C. Therefore, the frequency parameter depends 
only on frequency. No justification is given for evaluating k^ at a 
temperature of 0®C. The experimental results indicate that the 
radiation resistance and the reflection coefficient are not signifi- 
cantly dependent upon temperature for a given value of the frequency 
parameter; however, the resistance depends on temperature for a 
given value of the wavelength parameter. Agreement between experi- 
mental results and the theory of Ando^® is good when the radiation 
resistance and reflection coefficient are plotted against the fre- 
quency parameter. 

Cummings^® also uses the impedance tube -technique to investigate 
the effects of hot flow on the radiation impedance of an unflanged 
circular duct. Cummings' experimental results are similar to those 
of Pricker and Roberts in that they indicate a significant effect 
of temperature on the radiation resistance for a given wavelength 
parameter while showing no obvious trend in the radiation reactance. 
Cununings speculates that the normalized reactance represents a plug 
of air mainly at a temperature near that of the duct exit which 
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oscillates with the air in the end of the duct. Ho speculates fur- 
ther that the radiation resistance increases with temperature because 
the hot jot at the duct exit is equivalent to a continuation of the 
duct, with the temperature gradient between the duct exit and the 
ambient temperature acting as a horn to match the terminating im- 
pedance to the characteristic impedance of the surrounding medium. 

The results of Fricker and Roberts and Cummings have to be 
used with discretion due to the absence in both cases of an adequate 
description of the acoustical environment surrounding the open end 
of the duct. 

Cline^® uses the impedance tube technique to investigate the 
effects of hot flow on the radiation impedance of an unflanged cir- 
cular tube. Experimental results indicate that the use of the 
classical method to measure the acoustic radiation impedance of a 
tube with nonisentropic flow is in principle invalid because of the 
distortion of the standing wave pattern by axial temperature gradi- 
ents. However, the method is valid from a practical standpoint if 
the temperature gradients are not too extreme. Cline concludes 
that a valid method has yet to be demonstrated in this case, although 
he suggests an iterative technique in which an appropriate wave 
equation is solved subject to an assumed terminating impedance until 
good agreement is obtained between the computed and obsexrved stand- 
ing wave structure. 

Zorumski^® extends an equation derived by Morse for the radia- 
tion impedance of a vibrating rigid piston to give the generalized 
radiation impedances of all modes in circular and annular ducts 
which have arbitrary wall admittance and which terminate in an in- 
finite baffle. He presents an exact equation for u. generalized 
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impedances of axisymmetric ducts with an infinite baffle. The duct- 
wall admittance is shown to bo a significant parameter in the effect 
of termination reflections. 

Johnston and Ogimoto^^'^^ develop an approximation for the 
terminating impedance of a rigid, finite length, unflanged circular 
duct with negligible wall thickness. Uniform air flow inside and 
outside of the duct is accounted for in this approximation. Two 
main differences in the sound field between the finite and semi- 
infinite duct are noted. First, the far-field directivity pattern 
changes from a distorted spherical pattern to a lobed pattern. 
Secondly, the standing wave pattern inside the duct changes due to 
the added reflections from both ends of the pipe. 

Kinsler and Prey^^ discuss the transfer impedance of both a 
travelling v;ave and a standing wave. The transfer impedance of a 
free travelling plane wave along a duct is equal to the character- 
istic impedance of the medium divided by the cross-sectional area 
of the duct. In the standing wave case where both incident and 
reflected waves are piesent, the acoustic impedance varies from 
point to point due to the varying phase relationship between the 
two waves. Thus, the transfer impedance depends upon the termi- 
nating impedance, the length of the pipe and the wave number. 
However, the terminating impedance depends only upon the physical 
nature of the termination. 

It is apparent that noise generation from a combustion process 
is not well understood. However, some approximate solutions of com- 
bustion noise have been postulated. More experimental data and 
research is needed in this area in order to understand the source 
mechanisms of combustion noise. 


13 


3. DEFINITION OF THE ACOUSTIC SOORCE/PROPAGATION MODEL 


The acoustic source/propagation model is based oh a first-order 
perturbation analysis of the one-dimensional unsteady flow of an 
inviscid ideal gas. The gas is assumed to be thermally nonconduct- 
ing, although heat loss to the burner walls is allowed through use 
of a negative volumetric source term in the corresponding steady 
flow energy equation. There is no attempt in the model to balance 
chemical species; combustion is treated as a locally known volnmetvic 
heating rate in which the chemical composition of the flow remains 
that of theoretical air, 

•j" 

Because the equations which describe the acoustic field in the 
burner are linear, the dependent variables can be related at any fre- 
quency of interest by a system of linear algebraic equations; 


■'’ll' 


’“ll “l 2 ^13 “14' 

^II 

a 

“21 “22 *^23 “24 

r 

H 


^31 °'32 “33 “34 

km tm 


Q 


( 1 ) 


where the indices I and II refer to any two sections of the combustor 
where the complex values of pressure P, velocity V and temperature T 
are defined. 

The quantity Q in equations (1) is a complex number whose value 
depends on frequency and which represents the ratio of the unsteady 


We refer to the unsteady component of the flow field in the burner, 
represented mathematically by the perturbation variables , as an 
"acoustic” field, even in the source region. We recognize that 
such a designation may be pedagogically incorrect, but it is con- 
sistent with the idea of an "acoustic” model. 

This "transfer matrix" approach was suggested by E. A. Krejsa^® on 
the basis of a similar analysis by J. H. Miles^°. 
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heat release at that frequency to the total steady heat release. 
Thus it is the thermal-acoustic efficiency whose distribution with 
frequency we refer to as the acoustic source spectrum. In repre- 
senting the source term by a single complex number at a given 
frequency we are assuming that the local acoustic sources are pro- 
portional to the local J.teady heat release at each frequency. ThJ,s 
hypothesis is based on a series of reasonable assumptions, the 
validities of which remain to b<( established; 

1. Only the local fluctuations of h<6at release are dom- 
inant in producing the local fluctuations of density 
which give rise to the acoustic waves propagating 
into the far field as sound (l.e., we consider only 
the possibility of direct combustion noise) . 

2. The only significant fluctuations of heat release 
are those due to the local turbulent mixing of the 
unburned fuel and air with the hot products of com- 
bustion. 

3. The local turbulence level (integrated over the 
burner cross-section) is uniform in the combustion 
zone. 

Convincing arguments in favor of the first two assumptions are 
presented by Legendre in reference 22 and will not be reiterated 
here . Although the third assumption is somewhat more questionable , 
this does not impose a serious limitation on the results obtained 
if a realistic form is assumed for the axial distribution of steady 
combustion heat release. For example, in the present study we 
assume a distribution of the form 


q(x) ” 1 - cos (21TX/L)] , 


(2) 


where q^ is the maximum value of the local steady heat release and 


L is the length of the combustion zone. This distribution is similar 
to that actually found in the experimental burner. Moreover, Mahan 
and Kasper 3 6 have established that the low frequency acoustic response 
of a long burner is not very sensitive to the heat release distribution 


as long as the oorrect total amount of heat is released in the com- 
bustion zone, Then for the third assumption to be valid from a 
practical standpoint it is only necessary for the turbulence level 
to be uniform in the relatively limited region of the combustion 
zone surrounding the peak in the distribution given by equation (2) 
in which most of the heat is released. 

The elements a. , of the transfer matrix of aquations (1) are 
ij 

determined by solving the equations which describe the unsteady one- 
dimensional flow of an inviscid ideal gas with a distributed heat 
source : 

the continuity equation, 
the equation of motion, 

pv^ + pw^ + p^ « 0 , (4) 

the energy aquation, 

pCpTt + “ Pfc 

the thermodynamic equation of state , 

p w pRT , (6) 

and the calorimetric equation of state, 

c - c (T) , (7) 

P P 

where the variation of specific heat with temperature is computed 
assuming a mixture of 79 per cent nitrogen and 21 per cent oxygen 
and using the low pressure formulas given by VanWylen and Sonntag^^. 
If one makes the usual assumptions that the solutions p, v, T and p 
consist of a steady component which is a function of axial position 
only and a very much smaller fluctuating component which is a func- 
tion of both axial position and time, and that the steady solution 
itself also satisfies equations (3) through (7), there result 
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v'/v + p'/P “ T'/T •• 0 , 

pvv* + RTp' " 0 

and pvCpT' “ RTvp' + R^'q ; 

and 

(e 5 - T)^ + ^(p + \? " T)^ « 0 , 

<}. + (2^> - T)'7' + + (RT/^p » 0 

vi X X 

and + v*T^ « [p^ + - (p + ^ - q)q/p] . 

In equations (8) and (9) - the quantities with overbars ("~) are the 
steady flow variables and the quantities with tildes ("') are the 
acoustic variables nondimensionalized by the corresponding steady 
flow variables. The subscripts t and x indicate partial differen- 
tiation with respect to those variables and the prime ( ' ) associated 
with tho steady flow variables indicates differentiation with respect 
to axial position. 

Equations (8) can be integrated directly for a specified axial 
distribution of steady heat release q(x) beginning at any section, 
the burner inlet for example, where initial values of the pressure, 
temperature and velocity are known. Once these equations have been 
solved the varieible coefficients in the acoustic equations are known. 
In order to solve the acoustic equations (9) , periodic solutions are 
assumed of the form 



iwt 

^ = Pe 
■* „ iwt 

and 

m iwt 

•p S3 *p0 

with 

„ _ io)t 

§ = Qe 


where the quantities P, V, T and Q are complex functions of axial 
position only. When these assumed forr«(S are introduced into equations 
(9) there result the ordinary differential equations 
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iw(P - V) + V(P' + V - T') - 0 , 

iwV + (2V - T)v' + W + rTTp'/v * 0 


and 


iwT + vlf* ■ (*^*^) ’*■ "* ^ Q)¥/p] • 


ill) 


which can be integrated numerically from any section where initial 

I 

t 


values are known. 

Unfortunately,, a straightforward integration of the acoustic 
equations (11) is not possible because, unlike the case of the steady 
flow equations (8) for which q is specified and the values of p, v 
and T are known at the burner inlet, there is no section where the 
complex values of P, V and T are known and, of course, it is the 
variation with frequency of the thermal-acoustic efficiency Q that 
is sought in this study. It is for this reason that we have adopted 
the transfer matrix approach represented by equations (1) . We can 


find the elements of the matrix by integrating equations (11) 


four times using each time a different but linearly independent set 
of arbitrary initial cdnditions P^, V^, and the source function 
Q. If, for example, for the first integration the initial values 


at the inlet of the burner are assumed to be » 1 , = 0 


and the source function is assumed to be Q « 0, the first column 
of coefficients can be expressed in terms of the solutions obtained 
by integrating along the burner to the open end of the exhaust pipe; 


i.e., = P„, = V,, and a 


Tjj. Once the elements of the 


II-' "21 'II 31 

transfer matrix have been thus deteimined for a given frequency. 


there remain three algebraic equations in seven unknowns : the 

values of P, V and T at each end of the burner end the value of Q. 


t . 


Since P, V, T and Q have real and imaginary parts, equations (11) 
really represent six equations. However, for purposes of the 
discussion which follows we treat this as a system of three equations, 


rnimaairriiiriffiimiim 




MmMSifiiffiiriiaim 
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Three additional equations can be written in terms of the measured 
acoustic impedance at either end of the burner, 

Zj - Pj/Vj and Zjj - , (12) 

and the isentropic condition at the burner inlet, 



The final equation needed to have a determinant system is the relation 
between the pressure fluctuation at the open end of the exhaust pipe 
and the acoustic pressure measured at the same frequency in the far 
field; 

P__ = A{ai)P^ , (14) 

XX o 

where A(u) is the appropriate transfer function. In principle, we 

are now in a position of being able to compute Q using equations 

(1) , (12) , (13) and (14) for any frequency co for which we have 

measured P . 

o 

The foregoing analysis pertains to the simple burner configu- 
ration without secondary air flow or holes in the liner, A brief 
discussion of the procedure which must be followed in the case of 
holes and secondary flows may be found in subsection 6.4. Discussion 


of Recommendations . 
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4 . IMPLEMENTATION 


The experimental apparatus is shown in Pig. 1. The fuel^ 
hydrogen, enters the combustion chamber through a small sonic ori- 
fice in the center of a disk around which air is introduced with 
a swirl component of velocity. The use of a sonic orifice assures 
that the acoustic-scale pressure fluctuations in the burner cannot 
modulate the injection of hydrogen, thus avoiding the possibility 
of combustion "feed” instability. The air, which flows through a 
silencer before entering the combustion chamber, is metered by a 
calibrated rotameter. The wake created by the disk and the swirl 
component of , velocity stabilizes the flame. 

The products of combustion are vented into an anechoic chamber 


through an exhaust pipe having the same inside diameter, 54.6 ram, 
as the combustion chamber.- The burner, including the exhaust 


pipe, is 1.53 m long and penetrates 0.78 m into the anechoic 
chamber. The flanges which couple the fuel and air injection mani- 


fold and the exhaust pipe to the combustion chamber are water 
cooled, as is that poi.‘tion of the exhaust pipe outside the anechoic 
chamber. The anechoic chamber measures 2.3 u x 2.6 m x 4.0 m and 
is ventilated by a silenced blower whose capacity limits the power 
level of the burner to a maximum of around 20 kW. The burner power 
level is computed from the hydrogen flew rate measured using a 
calibrated rotameter and assuming complete combustion. 

Two Bruel and Kjaer 4133 1/2-in. (1.27 cm) microphones are 
located in the anechoic chamber at a distance of 1.5 m from the 


open end of the exhaust pipe, one on the axis of the burner and 
the other at an angle of 45 deg in the vertical plane above the 
axis. These microphones are connected to B and K 2603 amplifiers 
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through B and K 2801 preamplifiers. A Zonic 6080 Fast Fourier Trans- 
form analyzer interacts with a Tektronix 4052 graphics terminal to 
produce autospectra from the sound pressure signals measured by the 
microphones. The system is calibrated using a B and K 4220 piston- 
phone . 

4.1 Measurement of the Burner Terminating Impedances 

The classical impedance tube method was used to measure the 
acoustic impedance at each end of the burner. The measurements were 
made in the cold burner without flow and the results extended to the 
cases with combustion using the relation 

Z“pczv/p , (15) 

where pc is the characteristic impedance corresponding to the local 
conditions with combustion and z is the specific impedance measured 
without combustion. There is no need in this study to correct the 
impedance for flow effects because the Mach number is always less 
than 0.02 in the burner. 

There can be little doubt of the validity of equation (15) at 
the burner inlet because the conditions there never differ greatly 
from those without combustion. However, the conditions at the out- 
let of the exhaust pipe Vciry significantly with power level. Indeed, 
the rather steep axial temperature gradient which exists at the end 
of the burner must cast doubt on the validity of using this method 
to determine the impedance there. It is possible, for example, that 
these gradients can alter the degree to which a circular pipe radiates 
acoustic energy. Two recent studies, discussed below, seem to indi- 
cate that this is not the case, although more research is needed to 


confirm this conclusion. 


I 
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Pricker and Roberts^^ and Cummings^® both uaed the claasical 
impedance tube technique to measure the radiation impedance of a 
thick-walled tube from which a hot gas flowed into a relatively 
cold space. In both experiments the flow velocity was too low to 
significantly influence the radiation J.mpedance. Both the studies 
show that although the imaginary component of impedance does not 
seem to be affected by the difference between the gas temperature 
in the pipe and that in the room, the real component, when expressed 
as a function of wavelength in the pipe, increases with the tempera- 
ture of the gas in the pipe. However, when the real component is 
expressed as a function of frequency rather than wavelength it is 
relatively independent of the gas temperature in the pipe. This is 
exactly the behavior predicted by the classical isothermal theory of 
Ando^® for a semi-infinite pipe opening into an anechoic space.. 

The work of. Pricker and Roberts and Cummings then seems to indi- 
cate that equation (15) is valid for characterizing the radiation 
impedance at the open end of the exhaust pipe for a given frequency 
if it can be assumed that the impedance tube method itself is appli- . 
cable in the case of axial temperature gradients. That this is the 
case is certainly open to question because the standing wave structure 
in the tube, assumed to be uniform in the analysis associated with 
the classical impedance tube technique, is greatly distorted by the 
presence of axial temperature gradients. In order to study this 
question we have developed an impedance measurement technique which 
is independent of the axial temperature distribution in the exhaust 
duct. We have utilized this technique in a preliminary study to 
verify the measurements obtained using the classical method. 
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a?he apparatus for this preliminary study is shown in Fig. 2. 
It consists of a Pyrox tube 1,0 ro long and having inside and out- 
side diameters of 74,6 and 79,8 mm, respectively. The lower end 
of the tube is connected to a small settling chamber and the 
upper end opens into the laboratory. A spenker located in the 
settling chamber and driven by an oscillator produces a standing 
wave in the tube at a specified freg;uency. Near the lower end 
there is a small natural gas burner which can be sUsplaced along 
the axis of the tube. Air enters the tube from the settling 
chamber and mixes with the natural gas, the mixture is burned, 
and the products of combustion exit from the open end. 

The acoustic pressure in the tube downstream of the burner 
is measured by a probe borrowed from a B and K 4002 impedance 
tube apparatus. The probe must bo introduced from the upper end 
of the tube because the burner occupies the axis on the other 
end. Great care has to be taken to assure that the temperature 
of the probe is constant during the measurements and that the 
temperature is the same for each measurement. If not, the sensi- 
tivity of the probe will be different from one measurement to the 
next, thus introducing significant scatter in the data. We were 
able to avoid this potential source of error by plunging the 
probe rapidly into the tube and obtaining the pressure measure- 
ment before the temperature of the probe could change. After 
each measurement the probe was retracted and allowed to attain 
its original temperature before the next measurement. 

In addition to the local ’acoustic pressure on the axis of 
the tube, the temperature of the gas at the exit of the tube and 
in the laboratory as well as the air and natural gas flows were 


measured. 


The usual ajpproach in the standing wave, or impedance tube/ 
method consists of measuring the maximum and minimum values of 
the acoustic pressure amplitude in the tube as well as the axial 
position of the pressure node nearest the open end. Then the 
radiation impedance is computed from the standing wave ratio and 
first node position under the assumption that the standing wave 
pattern is uniform throughout the tube. In general this is true 
only if the temperature of the gas in the tube is uniform; if there 
is a significant axial temperature gradient/ the values of the 
maxima and minima as wall as the wavelength of the standing wave 
will vary along the tube/ thus jeopardizing the validity of the 
method. Pricker and Roberts tried to avoid this problem by 
extrapolating the standing wave ratio measured at several locations 
along the tube to the end of the tube/ but they evidently did not 
perform a similar correction for the wavelength. Cummings / on 
the other hand, tried to avoid axial temperature gradients alto- 
gether by selectively heating the walls of the tube. 

Our approach consists of finding the radiation impedance/ 
for example by iteration/ which minimizes the r.m.s. difference 
between the measured and calculated standing v/ave patterns. The 
calculated pattern is obtained by first integrating equations (8 ) , 
which describe the mean floW/ assuming a uniform heat loss to the 
walls chosen such that the calculated gas temperature at the tube 
exit is the same as the measured value. Then the acoustic equa- 
tions (11) are integrated starting at the open end of the tube 
with the measured value of the acoustic pressure p| with a value 


+ 

This value was slightly adjusted in several cases to improve the 
agreement between the measured and calculated waveforms / thus 
avoiding according too much importance to a single member of a 


t, 
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of the particle velocity V given by equation (12) / and with a 
value of the acouetic temperature T given by equation (13) . 

Because the integration was always terminated before reaching 
the flame, the value of Q was taken to be zero. 

Figure 3 shows the variation of the real and imaginary parts 
of the radiation impedance as a function of the wave number 
based on the temperature of the gas in the tube and reduced by 
the internal radius of the tube for a range of gas temperatures 
at the tube exit. These results were obtained by applying the 
classical analysis, which assumes uniform (undistorted) standing 
waves, to the data obtained as described above. That is, the 
standing wave ratio and the axial position of the first node 
have been used without correction for the distortion of the 
standing wave by the axial temperature gradient. The impedances 
thus obtained are in qualitative agreement with those of Pricker 
and Roberts and of Cummings in that the real part increases with 
the temperature of the gas in the tube for a given reduced wave 
number. The significant differences which exist between the 
experiments of Fricker and Roberts and Cummings and our own 

+ 

apparatus preclude the possibility of a quantitative comparison. 

Figure 4 shows the measured and calculated pressure standing 
wave waveforms for five cases representing a range of heating con- 
ditions. In each case the frequency was 1000 Hz, the air entered 
the tube at a temperature of 17 ‘’C and with a mean velocity of 
0.25 m/s, and the temperature in the laboratory v/as 17 ®C. The 
thermal conditions in the tube were varied by changing the position 

i* 

Indeed, there were significant quantitative differences between 
Pricker and Roberts ' results and those of Cummings . 


I 


as 


and boating lavol of the burner and, in one caio, by wrapping the 
tube with thermal insulation. The symbols in the figure represent 
the measured values and the curves are the theoretical waveforms. 
The theoretical waveforms wore obtained by integrating the acoustic 

equations starting with initial conditions based on impedances 

* 

taken from Fig* 3. The good agreement between the experimental and 
theoretical results in Fig, 4 demonstrates that good values of 
impedance can be obtained by analyzing the data in the classical 
way without considering the distortion of the wave. Thus we can 
conclude that the classical impedance tube method used by Fricker 
and Roberts and by Oommings, if the data are obtained carefully, 
gives good results even in the presence of axial temperature 
gradients, at least for the limited range of frequencies and 
thermal conditions of this study. Then it seems that we can in- 
deed accept their results and use with confidence equation (15) 
for converting the impedance measured without /tjnmbustion to an 
impedance corresponding to combustion conditions at a given fre- 
quency . 

The impedances measured at the two ends of the burner are 
shown in Pigs. 5 and 6. The theoretical impedance of Ando also 
appears in Fig. 6 , which shows the magnitude and phase of the 
specific impedance at the open end of the exhaust pipe. The agree- 
ment between experiment and theory is very good above 700 Hz. The 
divergence which begins below 700 Hz is probably due to the depar'** 
ture of the chamber from anechoic conditions at low frequencies. 

4.2 Determination of the Source Spectra 

The strategy for solving equations (1) with equations (12) 
through (14) is as follows. First, the steady flow equations are 
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integrated subject to known initial conditions and a specified axial 
heat transfer distribution. In principle these include the pressure, 
velocity and temperature at the inlet of the burner, the known com- 
bustion heat release distributed according to equation (2) , the 
measured heat loss to the flanges and the cooled portion of the ex- 
haust pipe, and the value of the heat loss to the uncooled portion 

i 

of the exhaust pipe necessary to obtain the measured gas temperature 
at the open end. In practice, however, a better agreement can be 
obtained between the frequencies at which peaks occur in the cal- 
culated and measured press ..e spectra by slightly adjusting the 
assumed heat loss to the uncooled portion of the exhaust pipe. In 
fact, the distribution of the pressure peaks with frequency is 
fairly sensitive to the distribution of heat loss along the burner 
for a given total value. An obvious improvement to the technique 
would be to measure this distribution very carefully and with better 
spatial resolution. 

Next an arbitrary value for Q is assumed whose magnitude is 
constant at The six equations in six unknowns which remain after 

the value of Q has been specified are then solved for each frequency 
of interest. For these calculations the flame is divided axially 
into 20 source elements. Then in order to have incoherent sources, 

20 spectra were calculated corresponding to a sowce at each of the 
source locations along the length of the flame. These 20 spectra 
were then added together to form what we call the primitive spectrum. 
The use of incoherent sources is consistent with the idea that the 
sources are due to turbulence. 

After having thus obtained a primitive pressure autospectrum 
whose peaks align fairly well with those of the measured pressure 
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autoapdctjrum, tha corraapondlng aourca apactrum can ba calculated 
from 

Q(w) « Pjj(w)Q^/P(w) , ( 16 ) 

what. 1. obClnad from th. m...ur.d aound .p«otr« u,ln, 

equation (14) and P(m) la the calculated primitive aound apectrum 
corraapondlnq to Q^. The validity of equation (16) la baaed on the 
fact that equationa (1) are linear. In the reaulta preaentad In 
this report a value of unity haa bean aaaumod for the tranafar 
function A(w) in equation (14), Although this function was not 
measured, the transfer function between a location 20 cm inside the 
Qxhaust pipe and the microphone location was measured. This func- 
tion was constant at a value of about 0.75 over the frequency range 
from 150 Hs to 600 Hs, beyond which a meaningful measure could not 
be obtained because of the presence of a pressure node near the 
microphone inside the pipe, A value of unity for A(w) seams justi- 
fied in view of the relative proximity of the field microphone to 

♦ 

the open and of the pipe and the limited frequency range. 
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5 . RESULTS 


S . 1 The Measured Specfcra 

Typical background and air flow noise spectra without combustion 
corresponding to the three combustion liners used (see Fig. 7) are 
shown in Pigs. 8, 9 and 10. Liner I has no holes, liner II has a 
ring of four evenly spaced 12.7 mm diameter holes 76.2 mm from the 
inlet end of the burner, and liner III has a ring of four evenly 
spaced 12.7 diameter holes 228.8 mm from the inlet end of the burner. 

Figures 11 through 16 compare the spectra at the two microphone 
locations for three different power levels, 4.5, 10.0 and 22.3 kW, 
for each of the three liners. In all cases in these figures the 
total air flow is 620 standard liters per minute (SLPM) , and in the 
cases of liners II and III the primary flow through the swirler is 
470 SLPM and the secondary flow through the holes in the liner is 
150 SLPM. As expected, the 22.3 kW power level always produces the 
highest spectrum, followed by the 10 kW power level and finally by 
the 4.5 kW power level. However, the dependence of the acoustic 
pressure on thermal power is not linear for any of the combustion 
liners, indicating that the overall thermal-acoustic efficiency of 
the combustion process is not constant. The noise spectra for the 
three different power levels converge at high frequencies (~1500 Hz) 
to the level without combustion (compare with Figs. 8, 9 and 10) for 
all three combustion liners, indicating that the effects of combustion 
are limited to low frequencies. Of course, the peaks in the sound 
spectra correspond to the longitudinal modes of the burner. These 
modes are rather diffuse for combustion liner II; in fact, it is 
easy to imagine that there are double peaks , especially at the lower 
frequencies. This is probably due to the steep temperature gradients 
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associated with the introduction of secondary air in the flame zone. 
This possibility is discussed later in this section in connection 
with the explanation of double peaks observed in the theoretical 
spectra. The distance between the peaks increases as the power 
level increases. This is of course due to the temperature depend- 
ence of the speed of sound. Comparisons between the sound spectra 
at the two microphone locations show only minor differences at all 
power levels and for ill three combustion liners. Thus we will 
limit our study to the spectra obtained from the "on-axis" micro- 
phone . 

Figure 17 compares the effects of two different air flow rates, 

290 and 620 SLPM, for combustion liner I. The longitudinal modes 
for this combustion liner are very distinct. The higher air flow 
rate case shows much higher peaks , suggesting that the increased 
level of turbulence associated with the higher Reynolds number flow 
enhances the thermal-acoustic efficiency. The lower air flow rate 
case corresponds to higher mean flow temperatures along the burner 
and exhaust pipe , thus producing larger distances between the 
resonant peaks. 

Figure 13 compares the effects of three different air flow com- 
binations for combustion liner II. As pointed out in connection 
with Figs. 13 and 14, the peaks corresponding to the longitudinal 
modes are not as distinct with this liner as with combustion liner I. 

The case for which all of the air enters through the swirler has the 
highest sound pressure level. This is perhaps because the flow Reynolds 
number at the base of the flame is largest in this case. Even though 
another case has the same total air flow, the sound pressure level is 
lower because part of the air flow is introduced through holes in the 
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liner wall, thus reducing the mean Reynolds number and presumably 
the turbulence level at the base of the flame. The "jet'' velocity 
through the holes in the liner is much too low to contribute to 
the turbulence of the flame. The case with the lowest total air 
flow rate produces the lowest sound pressure level, due probably 
to the decrease in turbulence. 

Figure 19 compares various air flow combinations for combustion 
liner III. In the frequency range between 500 and 700 Hz where the 
highest sound pressure levels are obtained, higher primary air flow 
rates lead to higher sound pressure levels. This lends weight to 
the argument that the dominant factor for influencing the thermal- 
acoustic efficiency is the Reynolds nximber in the flame region. As 
in all other cases we note that the spectra approach the "no-flame” 
levels of Figs. 8, 9 and 10 as we approach 1500 Hz. 

A direct comparison between the three combustion liners ,for a 
power level of 4.5 kW and a total air flow rate of 620 SLPM is shown 
in Fig. 20. Combustion liner I has the highest sound pressure level 
of the three liners , presumably because all of the air enters through 
the swirler. Combustion liner II and III have the same total air flow 
rate as liner I but a portion is introduced as secondary air through 
the holes in the liner. Comparing the sound pressure spectra for 
these latter two liners shows that each produces higher acoustic pres- 
sure levels in certain frequency ranges. The frequency interval be- 
tween resonant peaks is very nearly the same for each of the combustion 
chanibers , indicating that the steady-state temperatures are not signifi- 
cantly different. 

Figure 21 compares the sound pressure spectra produced by com- 
bustion liners I, II and III for a power level of 10.0 kW and a total 
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air flow rate of 620 SLPM. The results shown are very similar to 
those shown in Pig, 20, and thus the same observations can be made. 

The overall sound pressure level for all three liners has increased, 
as has the interval between resonant peaks, due to the increased 
power level and temperature, 

Figure 22 compares the spectra produced by combustion liners 1, 

II and 111 for a power level of 22.3 kW and a total air flow rate of 
290 SLPM. As noted previously in the discussion of the lower power 
level cases, combustion liner 1 produces the highest noise level over 
the whole frequency range due, once again, to the fact that all of 

the air is introduced through the swirler in the head end of the burner. 
Above 500 Hs the sound pressure level produced using combustion liner 

III is higher than that corresponding to combustion liner II* This may 
bo due to a "detuning" of the burner as suggested in connection with 
the discussion of Figs, 13 and 14. In this scenario one can imagine 
peaks which are broader but still contain as much or even more acoustic 
energy as the narrower, sharper peaks associated with combustion liners 
1 and 111, This is evidently the case, for exan\pla, for the peak around 
775 Ha in Pig. 22. 

Figure 23 compares the spectra produced by combustion liners 1, II 
and III for a power level of 22,3 kW and a total air flow rate of 620 
SLPM. Only minor differences occur between the spectra for combustion 
liners 1 and III for these, the highest levels of burner power and total 
air flow rate. However, the peaks as well as the valleys of the sound 
pressure spectrum corresponding to convbustion liner II are significantly 
lower. Secondary air introduced into the fl«une region (or the presence 
of holes there) clearly has the effect of detuning the burner rather 
than the anticipated affect of increasing the thermal-acoustic efficiency. 
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Figure 24 compares the spectra produced by combustion liners I, 

II and III for a power level of 10. 0 kW and with all of the 620 SLPM 

of air introduced through the swirler. The only parameter varied in 

this figure is the hole pattern in the combustion liners. The sound 
pressure spectra of combustion liner I dominates over most of the 

frequency range. This implies that, as might be expected, the pres- 

ence of holes in the liner tends to detune the burner in much the 
same way that an expansion chamber exhaust muffler attenuates sound 
from an internal combustion engine. The hole pattern of combustion 
liner II surprisingly has a less attenuating effect than that of 
combustion liner III. Thus comparison of Pigs. 23 and 24 leads us 
to conclude that air flow through the holes of liner II is an addi- 
tional factor for detuning the burner. 

5,2 Recovery of the Thermal-Acoustic Efficiency Spectra 

We now apply the data interpretation technique developed in 
Section 4 to the three sound pressure spectra shown in Fig. 11. The 
corresponding primitive spectra, calculated using a value of vT x lo ^ 
for the magnitude of Q^, are shown in Fig. 25, These calculations 
were limited to 150 at the low end of the spectrum because the termi- 
nating impedcmces at either end of the burner could not be measured 
with confidence below this frequency. 

There are four obvious differences between the measured spectra 
of Fig. 11 and the calculated spectra of Fig. 25. 

1 . The measured spectra generally decrease with frequency 
while the calculated spectra generally increase. 

2. The measured spectra have simple peaks while the cal- 
culated spectra have double peaks above a certain 
frequency. 

3. The calculated spectra have larger peak-to-valley 
differences (and sharper peaks) than the measured 
spectra. 
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4, The general di^pression in the spectra around 700-900 
Hz (depending on the burner power level) is much more 
pronounced in the calculated spectra than in the 
measured spectjra. 

The first difference noted above is a direct consequence of the 
difference between the actual source spectra in the burner and the 
assumed constant source spectrum whereas the other three differ- 
ences are artifacts of the acoustic model. 

The double peaks in the calculated spectrum occur because there 
are two acoustic lengths in the model which produce resonances in 
the frequency range of the calculations; the distance between the 
two ends of the burner and the distance between the open end of the 
burner and the steep temperature gradient associated with the com- 
bustion heat release. The sharp variation of the characteristic 
impedance pc associated with this temperature gradient results in 

significant reflection of acoustic energy at sufficiently short 
» 

wavelengths, thus producing secondary resonances. As noted pre- 
viously, this same effect could be produced by secondary flow into 
the burner through the holes in the liner for the case of liner II. 
Presumably a third family of peaks would appear at still higher 
frequencies above the resonant frequency of the short section of 
the burner between the closed end and the steep temperature gradi- 
ent, The fact that no obvious secondary peaks appear in the 
measured spectra indicates that the one-dimensionality assumed by 
the acoustic model is not strictly valid in the flame zone. Indeed, 
the flame does not fill the flame tube in the actual burner, thus 
leaving an annular region surrounding the flame through which 
acoustic waves can propagate without being influenced by the steep 
axial temperature gradient in the flame. 
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The calcuXated spectra have larger peak-to-valley differences 
(and sharper peaks) than the measured spectra because the only 
acoustic loss mechanism in the modal is the radiation resistance 
at either end of the burner. In the actual burner there are other 
loss mechanisms associated with viscosity and thermal conductivity 
which further limit the amplitudes of the resonant peaks. These 
secondary losses are much leas important in determining the values 
of the spectra at antiresonance, i.e., when the burner is detuned. 
Thus, the values of the source spectra calculated at antiresonance 
using equation (16) are much more representative than those cal- 
culated at or near resonance. 

The general depression, or dip, in the spectra especially 
evident at antiresonance around 700-900 Hz, depending on the burner 
power level, is a consequence of the Rayleigh criterion^®, the 
integral form of which requires that pressure disturbances of wave- 
length X be amplified according to the value of 

P(x) *q(x) dx , 

where p(x) is the local acoustic pressure variation and q(x) is the 
local unsteady heat transfer. The dot (•) indicates a scalar multi- 
plication of complex quantities; that is, the phase relation between 
p and q must be considered in evaluating the integral. A pressure 
disturbance of wavelength X will be amplified for positive values of 
the integral and damped for negative values while remaining unaffected 
for a value of zero. Then the depression in the sound pressure auto- 
spectra especially evident in the calculated spectra is due to the 
first alignment with increasing frequency of a standing wave pressure 
node with the peak in the heat release distribution, a situation which 
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minimizes the value of this integral. The fact that this effect is 
barely noticeable in the measured spectra indicates once again the 
inadequacy of a one-dimensional model in the flame zone; in reality 
only that portion of the acoustic wave in the center of the tube is 
subject to the amplification required by Rayleigh’s criterion. 

The source spectra calculated using equation (16) and the meas- 
ured and calculated sound pressure autospectra of Figs. 11 and 25 
appear in Fig. 26, The ordinate is in decibels referred to the 
burner power level; that is dB » 10 log^^CQ) where Q is the thermal- 
acoustic efficiency given by equation (16). The dashed lines repre- 
sent the actual source spectra calculated using equation (16) but 
ignoring the secondary peaks in the theoretical spectra, and the solid 
lines are smooth curves joining the values of the source spectra 
corresponding to the antiresonances of the sound pressure spectra. 
Following the above arguments concerning the differences between the 
measured and calculated sound spectra, the solid curves are more 
representative of the actual source spectra than the dashed curves. 
This is especially true at the higher frequencies where secondary 
loss mechanisms would be increasingly important. The "h\imps" which 
occur in all three spectra between 700 and 900 Hz are a result of 
the departure from one-dimensionality in the flame zone of the real 
burner as explained above and thus are not really significant. 

It is not surprising to see that the source spectra roll off 
rather steeply with frequency; this is consistent with the idea that 
the most important source mechanism is turbulent mixing of the un- 
burned fuel and air with the hot products of combustion. The fact 
that at 1500 Hz the source spectra are down nearly 20 dB from their 
values at 150 Hz justifies the decision to limit the study to this 
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fraqueucy range. All three source spectra decrease an order of 
magnitude between 150 Hz and the beginning of the influence of the 
hump around 600 Hz. This means that the moat significant source 
activity is limited to the same frequency range, 0 to 600 Hz, 
typically observed in jet and open flame noise studies. 

It is interesting to note that the three source spectra lie 
within a 5 dB band and in fact intersect at several frequencies. 

In general the low burner power case produces the lowest source 
spectrum, but the high power case does not produce the highest 
spectrum. That the low burner power spectrum rises above the high 
burner power spectrum beyond 1100 Hz is probably due to the onset 
of a second hump in the corresponding theoretical sound pressure 
spectrum for the low power case which occurs when the second pres- 
sure node is aligned with the peak in the heat release distribution. 
This second hump does not appear in the two higher power cases 
because the increased burner temperature associated with these 
cases pushes this hump out to higher frequencies. The combined 
uncertainty of the experimental and theoretical methods employed 
in this study is such that it is reasonable to conclude that the 
thermal-acoustic efficiency Q is relatively insensitive to burner 
power level. This efficiency has a peak value on the order of 10 
at 150 Hz and decreases about two orders of magnitude for a one 
order of magnitude increase in frequency. 

Another interesting aspect of these source spectra is that in 
all three cases there is a critical frequency near 500 Hz above 
which the resonant peaks in the corresponding sound pressure spectra 
produce sharp depressions in the source spectra (dashed lines) but 
below which the resonant peaks in the sound spectra translate into 
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sharp peaks in the source spectra. The sharp depressions have al- 
ready been explained in terms of the exaggerated resonant peaks in 
the theoretical sound pressure spectra, but how can the peaks in the 
source spectra below 500 Hz be explained? A possible explanation 
is that the mixing process is enhanced at low frequencies by the 
acoustic agitation of the flame. This is consistent with measure- 
ments of flame transfer functions reported in the literature^® 
which show that flames respond somewhat like amplified first-order 
systems, with a relatively flat response up to a critical frequency 
beyond which the response rolls off steeply with frequency. 
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6, CONCLUSIONS AND RECOMMENDATIONS 


6.i conclusions 


The following conclusions car, be tJrawn from the results of this 
stu4y ; 

1. Enhancement of the sound pressure level in the far- 
field of the experimental burner by the combustion 
process is limited tb frequencies below 1500 Hz. 

This is true for the combustion heat release range 
between 4.5 and 22.3 kW and for several hole patterns 
and air flow distributions in the burner. 

2. The sound pressure level in the "active" part of the 
spectrum (0-1500 Hz) increases with both burner power 
level and primary air flow rate. 

3. The presence of holes in the liner, with and without 
air flow through these holes, tends to attenuate the 
sound pressure level in the far field by comparison 
to the cases without holes, 

4. When there is no secondary flow through the holes in 
the liner, the attenuation of the sound pressure level 
in th/fi far field mentioned in 3. above is more pro- 
nounced for the liner whose hole pattern is located 
downstream of the flame zone. However, when there is 
flow through the holes the opposite is ture. 

5. The longitudinal mode frequencies are not modified by 
the inclusion of holes in the liner, except possibly 
for the case where the hole pattern coincided with 
the flame zone. In this case there is some evidence 
of secondary resonances associated with the acoustic 
length bounded by the hole location and the open end 
of the burner. 

6. The thermal-acoustic efficiency of the experimental 
burner is not very sensitive to the burner power 
level, especially at the higher power levels. 

7. The thermal-acoustic efficiency has a peak value of 
about 10“*^ at low frequencies ("*150 Hz) and decreases 
at an average rate of 20 dB per decade between 150 
and 1500 Hz. 

8. There is evidence that the turbulent mixing of the 
unburned fuel and air with the hot products of com- 
bustion is enhanced near resonancf; by acoustic agi- 
tation below a critical frequency around 500 Hz. 
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6.2 Dlscuaoion of Conclugions 

Conclusions 1 and 2 above tend to reaffirm the importance of 
turbulence in the flame zona for the conversion of thermal energy 
to acoustic energy. That is, enhancement of the far-field sound 
pressure field by combustion is limited to the range of frequencies 
normally associated with turbulence and jet noise. The rather steep 
roll-off of the sound spectrum, as well as the source spectrum, 
between 150 and 1500 Hz is consistent both with the roll-off of 
turbulent intensity observed in other studies and the decrease with 
frequency of the flame mobility reported in reference 39 and alluded 
to in conclusion 8. 

Conclusions 3 , 4 and 5 tend to indicate that the presence of 
holes opening into the annular plenum surrounding the flame tube 
produces a net attenuation not unlike that associated with an ex- 
pansion chamber exhaust muffler. The fact that flow into the 
flame zone through the secondary holes tends to attenuate rather 
than enhance the far-field sound level indicates that the detuning 
effect is more important than the increased possibility for tur- 
bulence associated with this secondary flow. In fact, the results 
seem to imply that the highest level of turbulence in the flame 
zone is obtained when all of the air is introduced through the 
swirler . 

The conclusions pertaining to the thermal-acoustic efficiency, 
numbers 6, 7 and 8, pertain only to the case of the burner without 
holes in the flame tube; the efficiencies for the cases with holes 
have not yet been calculated. 
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6 . 3 Recommendations 


The following recommendations are made on the basis of this 
study : 

1 * The analytical acoustic sourco/propagation model 
should be upgraded to include the losses associated 
with viscosity and thermal conductivity, 

2. The model should be modified to in some way account 
for the two-dimensionality of the flame zone, of 
course, an alternative would be to modify the 
perimental burner in such a way as to insure one- 
dimensional combustion. 

3. Calculation of the source spectra should be ex- 
tended to the cases of combustion liners with hole 
patterns and secondary air flows. 

4. The value of the transfer function A(w) in equation 
(14) should be carefully measured. 

5. In future experiments the gas temperature distri- 
bution along the burner should be measured and 
checked against the distribution calculated using 
equations (8) . The wall heat loss distribution 
q(x) should be tailored to insure good agreement 
between the measured and calculated axial gas tem- 
perature distributions . 

6. The techniques and ideas developed in this study 
should be extended to a gas turbine engine. 

6.4 Discussion of Recommendations 


The first three recommendations can be implemented without 

further experimental work. Indeed, we are presently in the process 

+ 

of implementing recommendation 3 at ONERA. We also have plans to 
implement the first two recommendations, but this will probably not 
be done until 1983-84. 

Extension of the analysis to the cases of the combustion liners 
with hole patterns and secondary flows, recommendation 3> requires 


1 * 

'ONERA = Office National d 'Etudes et de Recherches A^rospatiales , 
the French National Office of Aerospace Research where the Principal 
Investigator is spending the 1982-83 academic year on sabbatical 
leave . 
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only that oquatlonfi (1) bo and solvod for aaoh subsection 

of the burner bounded by the ends of the burner and the hole io- 
, cations, including the annular region surrounding the combustion 
liner. The general approach would otherwise be the same as for 
the cases without holes. The impedance of the holes themselves 
will be computed from classical theory as given for example in 
reference 40. 

Recommendation 4 can be implemented easily upon the return of 
the Principal Investigator to VPI in the Pall of 1983. However, 
it is expected that this transfer function will be very near to the 
assumed value of unity over the frequency range of this study. 

Recommendation 5 would require further experimental work with 
the same or a similar apparatus. This would only lead to a refine- 
ment of the present results without greatly changing the values 
already obtained for the source spectra. Thusc this recomniendation 
would only be implemented if a decision was made to continue the 
experimental program for other good reasons. 

Recommendation 6 should only be implemented after recommendations 
1, 2 and 3. In this case recommendation 4 concerning measurement of 
the transfer function A{w) would also be implemented, It is strongly 
urged that recommendation 6 be implemented, either at NASA Lewis or 
in the gas turbine research laboratory at VPI. 
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Pig. 1. Cut-Away View of Experimental Burner. 
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Fig. 2. Apparatus for the Preliminary Study. 
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Fig. 7. Hole Patterns of Combustion Liners I, II and III, 
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Typical Background and Air Plow Noise for 
Liner I at On-Axis Microphone. 
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Fig* 10. Typical Background and Air Flow Noise for 
Liner III at On-Axis Microphone. 
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Pig. 11. comparison of Sound Pressure Spectra for 
Combustion liner I, Three Power Levels, 
(On-Axis Microphone) . 
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Pig. 12. Comparison of Sound Pressure Spectra for 
Combustion Liner I, Three Power Levels, 
(Off-Axis Microphone) . 
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Fig\ 14. Comparison of Sound Pressure Spectra for 
Combustion Liner II, Three Power Levels, 
(Primary Air Flow =470 SLPM, Secondary 
Air Flow = 150 SLPM, Off-Axis Microphone) . 
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Pig. 15. Comparison of Sound Pressure Spectra for | 

Combustion Liner III, Three Power Levels, | 

(Primary Air Flow = 470 SLPM, Secondary I 

Air Plow = 150 SLPM, On-Axis Microphone) . | 
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Pig. 16. Comparison of Sound Pressure Spectra for 
Combustion Liner III, Three Power Levels, 
(Primary Air Flow = 470 SLPM, Secondary 
Air Plow = 150 SLPM, Off-Axis Microphone) 
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Comparison of Sound Pressure Spectra for 
Three Air Plow Combinations for Liner III. 
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Pig. 20. Comparison of Sound Pressure Spectra for 
the Three Liners (4.5 kW, Total Air 
Plow = 620 SLPM) . 
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Pig. 22. Comparison of Sound Pressure Spectra for 
the Three Liners (22.3 kW, Total Air Plow 
= 290 SLPM) . 
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Fig. 23. Comparison of Sound Pressure Spectra for 
the Three Liners (22.3 kW, Total Air Flow 
= 620 SLPM) . 
















Frequency, Hz 


Pig. 26. The Thermal-Acoustic Efficiency Spectra 
Corresponding to the Sound Pressure 
Spectra of Fig. 11. 





